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Abstract 

We derive a set of new formulas for various distributions in dilepton pro- 
duction via virtual photon bremsstrahlung from pseudoscalar mesons and 
unpolarized spin-one-half fermions. These formulas correspond to the leading 
and sub-leading terms in the Low-Burnett-Kroll expansion for real photon 
bremsstrahlung. The relation of our leading-term formulas to previous works 
is also shown. Existing formulas are examined in the light of Lorentz co- 
variance and gauge invariance. Numerical comparison is made in a simple 
example, where an "exact" formula and real photon data exist. The results 
reveal large discrepancies among different bremsstrahlung formulas. Of all 
the leading-term bremsstrahlung formulas, the one derived in this work agrees 
best with the exact formula. The issues of Mj'-scaling and event generators 
are also addressed. 
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I. INTRODUCTION 



Even if electromagnetic phenomena rank among the best understood in particle physics, 
their merging with hadronic processes often brings ambiguities that cannot be resolved on a 
purely theoretical basis. It is generally believed, however, that many of the reactions in which 
photons and dileptons are produced can be described as bremsstrahlung from incoming and 
outgoing charged particles, at least within a limited kinematic range. 

It is well known that the cross section for production of photons with very low energies 
is uniquely determined by the cross section of the corresponding nonradiative reaction . 
In the dilepton sector, the situation is less clear. Several different bremsstrahlung formulas 
have been proposed ||7|-[l4|, which, as will be demonstrated in this work, do not agree with 



one another very well. Some of these formulas fail to satisfy constraints implied by general 
principles such as Lorentz covariance and gauge invariance. Besides the formulas cited above, 
a few additional formulas exist that have been designed for specific kinematic regions |]T3,TB . 
They will not be considered here. 

The purpose of this work is to present a set of consistent formulas for various distribu- 
tions in dilepton bremsstrahlung from pseudoscalar mesons and unpolarized spin-one-half 
fermions. We consider the terms that are proportional to the square of the nonradiative ma- 
trix element (leading term approximation) and its derivatives (sub- leading approximation). 
We first derive the formula for the most general quantity, namely the double differential 
cross section in the momenta of leptons. The correct form for it has not been yet known 
even in the leading term approximation. Then we arrive at the cross section in dilepton 
mass and momentum. Its leading part differs only by higher-order terms in the dilepton 



four- momentum from one of the already existing formulas |T2 |. 

To stress the importance of using the correct virtual bremsstrahlung formalism, we com- 
pare our formulas to those already existing in the literature. The comparison is made first 
on general grounds and is followed by an application of the formulas to a simple physical 
process. All formulas are scrutinized from the point of view of Lorentz covariance and 
from what we will call a global variable test. 



The global variable test is based on the finding |TB|,[TP| that in a one-photon approxi 



mation, gauge invariance leads to the following relation for the inclusive^ differential cross 
section in global dilepton quantities M (dilepton mass) and q (dilepton momentum) 



r(M=) {^^\ (1.1) 

M 



dM^d^q ^ ^ \ d^q 



We will omit the word inclusive in what follows. All our cross sections for photon or dilepton 
production are semi-inclusive (integrated over the final state hadron momenta in a reaction with 
given number and types of final state particles). Some of the relations are more general and hold 
also for inclusive cross sections, which are given as sums of the semi-inclusive cross sections over 
all possible reactions with chosen initial particles [as, e.g., ( |1.1D , which is valid even for exclusive 
cross sections]. In some cases, the sum over all possible reactions should be supplemented to make 
the relation inclusive. 
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with the function T given by 




where a is the fine structure constant and /i the lepton mass. The rightmost quantity in 
Eq. (|LlD is called the cross section for virtual photon production. It does not have a direct 
physical meaning, as it is not experimentally accessible. Nevertheless, Eqs. (|1 . 1| ) and ( p..2|) 
can test the soundness of theoretical formulas, because they show the only two places in 
which the lepton mass /i may and must (unless neglected) appear. From a technical point 
of view, the virtual photon cross section is calculable more easily than the dilepton cross 
section |T^. Of course, if one needs the distribution in the momenta of leptons, a more 



involved approach is unavoidable. 

Furthermore, if an otherwise identical reaction exists in which a photon is produced 
instead of a dilepton (this need not always be the case, viz., Tr"*"-/?" e^e~\ the relation 




:i.3) 



must be fulfilled [0] with the differential cross section for real photon production on the 
right-hand side. 



Relations analogous to ( |1 . 1|) and (|1.3| ) can also be written for the corresponding quantities 



in other than two-body reactions, namely for decays and processes with more than two 



particles in the initial state |2y] 



In the next section, we derive the formulas appropriate for the production of very soft (low 
mass, low momentum) dileptons via virtual photon bremsstrahlung in reactions with charged 
pseudoscalar particles. Section |T| deals with the virtual bremsstrahlung from fermions. In 
Sec. |rV|, we show the additional approximations which lead to the various formulas that 
have appeared explicitly or as a part of more complex expressions in the literature. All the 
formulas are examined to see if they fulfill the Lorentz covariance and global variable tests. 
In Sec. 0, we introduce a simple theoretical model of the process (P —>■ 7T^7c~e~^e~ (which 
has not been experimentally investigated yet) based on a successful description of the 
recently observed p2| decay p° —>■ 7r+7r~7. The former process will be a testing ground for 



various virtual bremsstrahlung formulas in Sec. with the theoretical distribution from 
Sec. [V| serving as a reference. We summarize our main points and add a few comments in 
Sec. |VII| . Some related issues are discussed in the Appendices. In Appendix ^ we show how 
My-scaling transpires from the leading term virtual bremsstrahlung formalism. Appendix ^ 
deals with the "exact" formula for the p*^ — >■ 7r"'"7r~7* branching ratio. Appendix addresses 
the issue of photon and dilepton event generators that conserve energy and momentum. 



II. DILEPTONS FROM VIRTUAL BREMSSTRAHLUNG OFF PSEUDOSCALAR 

MESONS 

In this section, we assume that all charged particles are pseudoscalar mesons. The 
technically more involved case of virtual bremsstrahlung from spin-one-half fermions will be 
treated in the next section. 
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A. Leading term approximation 



Let us consider a 2 ^ n hadronic reaction 

a + 6^1 + 2H \-n (2.1) 

and denote its matrix element as Aio = M.o{Pa,Pb,Pi, ■ ■ ■ ,Pn)- Our aim is to find the matrix 
element Ai of the reaction 

a + b^l + 2 + --- + n + l+ + r, (2.2) 

in which a soft lepton pair with the four-momentum q = p_^.+p_ is produced in addition to n 
hadrons. The dominant contribution to the matrix element M. comes from diagrams where 
the virtual photon is attached to one of the external legs (see Fig. [^). The diagrams in which 
a virtual photon is radiated from internal lines give subleading contributions. This is caused 
by the nonvanishing virtuality contributions of the type (p^ — m^) to the denominators of 
newly emerging propagators. Such terms do not appear if one of the two particles attached 
to the electromagnetic vertex is real (p^ = m^). See below. 

Using the Feynman rules of pseudoscalar electrodynamics we can immediately write 
down the contribution to the matrix element Ai from radiation of an initial [x = a, b) 

= -eQ^Mo (2.3) 



and a final (z = 1, ■ ■ ■ , n) 



2pi-q + 



state particle. Above, e is the positive elementary charge, (Qi) is the charge of an initial 
(a final) particle, and M is the dilepton mass (g^ = M^). As is customary in this field, we 
have supposed that the nonradiative matrix element does not change when an incoming or 
outgoing momentum becomes "slightly" off-mass-shell; for example, 

Mo{pa - q,Pb,Pl, ■■■ ,Pn) = Mo{pa,Pb,Pl, ■ ■ ■ , Pn) = Mq . (2.5) 

The lepton part is given by 

= ^,u^'-\p~h,v^'-\p+) . (2.6) 

With this choice of L^, the matrix element for photon production is obtained by the sub- 
stitution ^ e^, where e is the photon polarization vector. Summing the contributions of 
individual diagrams, we arrive at 

M = eMoJ-L (2.7) 

with 
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2^ ^^O^ „ /\^2~'~Z^^« 



^^^^ 2p,-g-M2 ^-"'2p,-g + M2 
The four-vector J satisfies the important relation 



(2.8) 



(2.9) 



i=l 



which reflects charge conservation. Squaring the matrix element ( |2.7| ) and summing over 
the spins of leptons, we get 



The tensor Lny is defined by 



A straightforward calculation leads to 



(2.10) 



(2.11) 



(2.12) 



where we have introduced the four- vector / = p+ — p_ as the difference of leptons' momenta. 
Using the above relation we can rewrite (|2.10|) in the form 



E 1-^1' = IM 



2 327ra2 
' M2 



(2.13) 



Inserting this into the relation for the unpolarized cross section of the reaction ( ^.2| ) leads 
to 



da 



1 



4:EaEb\Va - VJ 



327ra^ \_j2 _J_ f, J^2 
M2 [ ^ "^V 



y,2E,(27r)3 2EK27r) 



(2.14) 



After neglecting the dilepton four-momentum in the argument of the 5-function and inte- 
grating over the momenta of final hadrons, we get 



a' 1 



d^p+(Pp^ Svr^ M2 



dan 



where 



dan 



1 



\M^\\2'nY5{pa + Pb-Y.P^)X{ 



d^Pi 



\ 2E,{2Tif 



(2.15) 



(2.16) 
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is the infinitesimal cross section of tlie reaction ( |2.1|) . Double differential cross section 
Ej_E_(f'a^'^^ /d?pj^d?p- is the most general quantity that characterizes the production of 
pairs of unpolarized unlike-sign leptons. Knowing it, we can find any other distribution (but 
not vice versa). 

A little exercise from relativistic kinematics provides us with the following general for- 
mula 



E 



a" 



JP + " d^p+d^p. 



(2.17) 



where E = E^ + is the dilepton energy and dCl^ is the solid angle element for positron 
momentum in the dilepton rest frame. Using ( p.l5| ) and (|2.17|) , we arrive at another formula 
of our virtual bremsstrahlung formalism 



E 



d' 



a 



a 



dM^d^q dn+ 327r4 



'1- 



4/i2 



dan 



(2.18) 



Finally, integrating over the positron angles in the dilepton rest frame, we obtain the differ- 
ential cross section in global dilepton variables 



E 



d^a 



a 



1 



dM'^d^q 127r3 



4/i2 



M2 



(2.19) 



This is just formula (|1 . 1|) with the virtual photon cross section given by 

'd^a^'\ 



E 



d^q 



M 



a 
4^ 



-J'] da, 



(2.20) 



Sometimes it is advantageous to express in terms of three-dimensional vectors. This can 
easily be done 0] using the relation Jq = E~^ J-q, which follows from Eq. (|2.9|) . We thus 
get 
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J^ = (Jxn)V-^(J-n)^ 



(2.21) 



where n = q/|q| is the unit vector in the dilepton momentum direction. 

To investigate the limit ( p..3| ) of (|2.20|) , it is sufficient to realize that due to (p.21|) , the 
part of that is proportional to q^ will not contribute in the limit M — > 0. What remains 
can be written as 



d^ 
d^q 



a 

4^ 



{-Jr) daa 



[2.22) 



with u; = |q| and 



, Pi 



CO Ej 1 — V,- ■ n 



(2.23) 
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To make the formula more compact, we have introduced the variable Q[, which is identical 
with the charge of final particles and acquires the opposite sign for initial particles. The 
unspecified sum runs over all (initial and final) hadrons and Vj = Pi/Ei is the velocity of 
the ith hadron. Eqs. (|2.22| ) and (|2.23|) combine to give the well known leading term formula 
for real photon bremsstrahlung p^j25[] . 

The central results of this subsection are the formulas (|2.15|) , (|2.18| ), and (|2.19|) for 
various differential cross sections of dilepton production via virtual photon bremsstrahlung 
in the leading term approximation. 



B. Next-to- leading term approximation 

When going beyond the next-to-leading order, the nonradiative matrix element is no 
longer considered to be immune against the changes in incoming and outgoing momenta. 
We write, instead of Eq. (|2.5|), 

Moip. - q) = Mo - q'^^ (2.24) 

if one of the incoming momenta changes, and 

-Mo(p. + g)=-Mo + g"^ (2.25) 

to account for a change in one of the outgoing momenta. In the above two equations, we 
suppress the momenta that keep their "nonradiative" values. If we now simply incorporated 
the "corrected" values (|2.24|) and ( p.25|) into the expressions (|2.3|) and ( p.4|) for the radiative 



matrix elements, and summed these up, we would obtain a non-gauge-invariant quantity. 
After replacing the four-vector by the virtual photon four-momentum, it would not van- 
ish. The reason is that we have not yet included the "contact" terms, which are generated 
from the strong interaction Lagrangian by the minimal electromagnetic interaction princi- 
ple. The latter says that the electromagnetic interaction terms appear as a result of the 
substitution 

p° ^ - eQg''^' , (2.26) 

where the index fi is to be contracted with the real photon polarization vector or the virtual 
photon propagator. To find the contact terms in our case, we make a formal expansion of 
■Moip'i — eQig°'^), where pi denotes any of the incoming or outgoing momenta. We thus find 
that the contact term associated with the ith hadron is 

= -eQ.f^ . (2.27) 

For the same reason as stated in the previous subsection, the radiation from internal lines 
will not contribute in this approximation either. Putting it all together, the radiative matrix 
element comes out to be 
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M = e K-L , 
where L is the four- vector defined by Eq. (^.61), and 



dPi 



Pi-Q 



9' 



(2.28) 



(2.29) 



We retained only the leading and next-to-leading terms in q. The sum runs over both 
incoming and outgoing hadrons and the four- vector J is given by Eq. ( p78| ) . To get the cross 
section for producing a pair of unpolarized leptons, we need the quantity 



J2 \Mf = 4:7ra H^^L 
with symmetric tensor L^^, defined by (|2.11|) , and 



Keeping only the g-terms of the same order as before, we easily obtain 



(2.30) 



(2.31) 



1 



= \Mo\ J^r + 



Q^Q'J d \M^ 



-Pi, 



X 



K 1 g'"^q^ - g'^^q'^ 



(2.32) 



To simplify the formula, we have again used the convention Q'j = —Qj for the incoming 
hadrons and Q'^ = Qj for the outgoing hadrons. The differential cross section in lepton pair 
momenta now reads 



a" 



a 



(Pp+d^p- Svr^ M2 



-J' 



M2 



[I-Jf] dao + 



Q^Q' 



M^{pi-q){p,j-q) 



Pi-qyi'-Pj 



I- Pi l-pj + M pi-pj] q 



da: 



(2.33) 



The notation is same as we met in Eq. ( p.l5| ), except for 



dao. 



d\M, 



.{27l)%Pa+P,-Y.Pk)X{ 



d^Pk 



(2.34) 



Using the same procedure as in subsection [11 A] , we arrive at the differential cross section 
in global dilepton variables in the form ( |1 . 1|) with the virtual photon cross section given by 



E 



'd^ 

. d^q 



M 



a 

4^ 



-J') dao + E 



Q^Q' 



iPi-q){pj-q) 



X 



{Pi-q)Pj - {Pi-Pj)q^] dao^tis 
We will comment on the zero photon mass limit of this equation in Sec. [m| . 



(2.35) 
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C. Decays 



The description of virtual bremsstrahlung in hadronic decays can be achieved following 
the lines sketched for the two-body reactions. We need only to change the number of 
incoming particles to one and replace the cross sections by decay widths. This follows from 
the similar structure of the relations between cross section, or decay width, on the one hand 
and the matrix element squared on the other. Another important factor is the universality 
(with respect to the numbers of incoming and outgoing particles) of Eqs. ( |2.7| ) and ( |2.8| ). 

Let us consider the decay 

a->l + 2 + --- + n + /+ + r. (2.36) 
For its differential decay width we can write, in the leading term approximation, 



d^T^^ 1 /• r ,o 1 



dTo . (2.37) 



The quantity 



,r„ = — lM,f i2.m,. - Ep.) n ^ (2.38) 



is the invariant decay width into an infinitesimal element of the momentum space for the 
decay 

a^l + 2 + --- + n. (2.39) 

The other leading-term formulas can be easily modified as well. The differential decay width 
in global dilepton quantities is 

If a nonradiative decay ( |2.39|) contains only two particles in the final state, we can proceed 
further to get 

nM-) ^-l^ I (-.P) ia,- , (2.41) 



To dW^dq* ' ' 47r2 E* 

where the asterisk refers to quantities in the rest frame of the parent particle and q* = |q*|. 
For a hypothetical decay into two particles and a virtual photon it means 
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III. VIRTUAL BREMSSTRAHLUNG FROM FERMIONS 



It has been shown in that for virtual bremsstrahlung from fermions, the terms pro- 



portional to in the expression ( |2.8| ) can be obtained in the same wayas in the real photon 
bremsstrahlung off fermions (see, e.g., |2^). In order to obtain the form of J analogous to 
( p.8|) , additional approximations are required. One has to discard some terms while keeping 
others of the same order in q. The justification of such a procedure is unclear. 

One may hope that if the summing (or averaging) over spins of hadrons is performed, 
the additional contributions will rearrange and the formulas identical to those pertaining to 
the pseudoscalar case will be restored. This would resemble a similar development in the 
next-to- leading terms in real photon bremsstrahlung [H. 

To investigate such a possibility, let as assume that in the nonradiative reaction ( p^.lD , a 
charged fermion-antifermion (e.g., proton-antiproton) pair is produced. To be more definite, 
we assign the index 1 to the fermion {Qi = 1) and 2 to the antifermion {Q2 = —1). Now, 
the matrix element exhibits the form 



(3.1) 



where T = T{pa,Pb,Pi, ■ ■ ■ ,Pn) is a matrix in the spinor space. To simplify notation, we will 
display its arguments only if they differ from the values just shown. Squaring the matrix 
element ( |3.1| ) and summing over the spin projections si and S2, we get 



l-Mc 



Tr 



+ m)r(p2 — m)r'] 



(3.2) 



where F' = 7°r^7'^. The quantity (|3.2|) determines the unpolarized cross section of the 
nonradiative reaction ( p.lD with the two spin-one-half hadrons in the final state. For a later 
use, let us notice that 

: Tr 



dpi 



r'7„r + T'ipi + m)— + —{pi+ m)T 



{P2 



m] 



(3.3) 



and 



d\M, 
dp2 



Tr < (pi + m) 



r7ar' 



dr 

dpi 



{p2 - m)r' + T{p2 



m] 



dV 
dp2 



(3.4) 



As a next step, let us consider the corresponding dilepton- producing reaction (|2.2|) . We 
will concentrate on virtual bremsstrahlung from fermions and, for simplicity, take all the 
mesons neutral. In addition, we neglect any anomalous electromagnetic interactions. The 
changes in the strong interaction matrix element will be incorporated by 



T{pi + q) = T + q' 



dr 

dpf 



(3.5) 



The contact electromagnetic interaction term associated with ith fermion line leaving the 
strong interaction core comes out as 
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— —eQi 



dT 

dp, 



(3.6) 



Using the Feynman rules for spinor electrodynamics, Dirac equation in momentum space, 



and the properties of the 7-matrices, we find the matrix element of reaction (^l2 



(3.7) 



The four- vector has the same meaning as before [see (|2.6| )]. We have, neglecting higher 
than linear g-terms in the numerators. 



1 



4pi-g 4:p2-q 
In accordance with ( p.8|) , we denoted 

(2pi + qY 



^ ^[7^ q] + ig^y^ - g^^q-) E Q.^l^ • (3.8) 



i=l,2 



Pi-qdp: 



(2p2 + qT 



(3.9) 



At this point we can clearly see the difference between pseudoscalar and fermion cases. 
When we neglected the changes in the strong interaction core as well as the contact terms 
in pseudoscalar case, we obtained immediately the leading term approximation in the form 
( p.7| ). For fermions, the extra terms (those with commutators) prevent us from reaching the 
same goal. 

The sum over the spins of the matrix element ( |3.7| ) squared assumes the form 



(3.10) 



The tensor L^^, is given by Eq. ( p.l2|) . Up to the leading and next-to-leading order in q, 

H^"" = \Mof J^r + ^{ J'^ Tr[(pi + m) F {p2 - m) K'"] 

Tr[(pi + m) K" (pa -m)V']+ fi ^ z/} (3.11) 

with K''^ = 7o(-ft''^)^7o. A straightforward manipulation guides us to an expression in parts 
of which we are able to identify the right-hand sides of Eqs. ( p.3| ) and ( |3.4| ). After replacing 
them by corresponding derivatives we get 



H 



Pi, 



d\M, 



X 



'^t^ {p^■q){p^qy dpf 

P^ {g^y^ - g'^q"") + P", {g'y^ - g'^q' 



(3.12) 



This is identical with what we would get from Eq. (|2.32|) for two outgoing charged mesons, 
with only one difference. A simple square of nonradiative mesonic matrix element is re- 
placed by the sum over fermion spins. As an independent check of our result, we explored 
also other situations (an incoming fermion-antifermion pair, one incoming-one outgoing 
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fermion/antifermion) and reached the same conclusion. For the initial state fermions, the 
sum is replaced by the average. The generalization to more than one fermion pair is obvious. 

The tensor H^^^ is the central object of all bremsstrahlung formulas. We have therefore 
proven that the dilepton production via virtual bremsstrahlung off unpolarized spin-one- 
half fermions is governed, in the leading and next-to-leading approximation, by the same 
formulas as that off pseudoscalar mesons. 

Especially, if we neglect the terms proportional to derivatives of the unpolarized nonra- 
diative matrix element squared, the tensor H^'^ reduces to 



with given by Eq. ( p.8[ ). The key leading term approximation relations ( p.l5| ), ( p. 181 ) 



and ( 2.19|) , which are of most practical interest, are valid also for virtual bremsstrahlung 



from unpolarized fermions. 

It is a good check that our expression ( |2.35| ) for the virtual photon cross section meets, 
in the limit of zero photon mass, the unpolarized photon cross section calculated from the 
Burnett and Kroll [Bl matrix element. 



IV. A SURVEY OF VIRTUAL BREMSSTRAHLUNG FORMULAS 

A. Riickl formula 

In his work 0, Riickl suggested the formula 

^+^-3n:- = — T7^(^^ , (4-1) 



which links the cross section for production of dileptons via virtual photon bremsstrahlung 
to the bremsstrahlung cross section for real photons. The meaning of the symbols is as 
follows: p+ and p_ are the momenta of leptons, E± = {p_^ + fi^Y^^ are their energies, 
= + p^y is the dilepton mass squared, q is the momentum of photon, u = |q| is 
its energy, and a is the fine-structure constant. Because of the vanishing photon mass, it is 
impossible to satisfy the relation u = Ej^ + E^ simultaneously with q = p+ -|- p_. 

Let us investigate now how Eq. ( [4. ID copes with general principles we mentioned in the 



Introduction. The principle of relativity |jT^ requires that any meaningful formula must 



be relativistically covariant. This implies that if we view the bremsstrahlung process from 
another (primed) frame, we must find the same relation among the transformed quantities 
as we did in the original frame: 

. (4.2) 



d^p'+d^p'_ 27r2 M2 \ d^q' 



q -P++P 



We have assumed for simplicity that the velocities of colliding particles are collinear and 
performed a boost along the collision axis (otherwise ex's also acquire primes). The left-hand 
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sides of Eqs. ( [4.1|) and ([4.2|) are obviously equal. Relativistic covariance will thus be satisfied 
if and only if 

But this condition cannot be fulfilled because in the new frame the photon momentum 
differs from the sum of lepton's momenta. In fact, for the longitudinal components of the 
corresponding vectors we have 

V±^L = 7 {P±,L - f3E±) 
q'L = l{qL-Pu;), (4.4) 

which leads to 

i'l = P'+,L + P'-,L + Pi (^+ + E_-u) (4.5) 

with a nonvanishing extra term on the right-hand side. 

In order to apply the global variable test, let us first use Eq. ( ^.17] ) to cast the Riickl 
formula in the form 



dM^d^qdQX 87T^M^\ d^q ' ^ ' ^ 

Integration over the positron momentum angles is simple because nothing depends on them: 



d'^a" " _ a 1 /. _ V d^(^^ 

dMH^q ~ 2^M2 V 'JP^^l^' ^ ^ 

This formula does not obviously have the form required by Eq. ( |1 . 1| ) . If we nevertheless 
extract from it the virtual photon cross section, which is defined by ( |1 . 1|) , we arrive at 

\W I M ~ 2 + M2 M2 + 2/i2 d^q • 

The M ^ limit of this expression is zero. If we are not so strict and require only 
/X -C M ^ C(j (this is the situation met, e.g., in the low-mass, high-transverse-momentum 
dielectron production), we are left with another surprising relation 

(4 9) 

The presence of an incorrect numerical factor in the Riickl formalism was already noticed 
by Craigie [T^, who ascribed it to the unjustified omission of the term Z^Zi, in the lepton 
tensor [see Eq. (|2.12|) ]. 

It is clear that the principal problems we have discussed here are not germane to the Riickl 
bremsstrahlung formalism, but are common to all the approaches where the dilepton cross 
section is assumed to be proportional to the photon cross section at the same momentum. 

To obtain Eq. ( [4.1| ) in our formalism, we must (i) write Pi-q = |q|-Ei — p.j-q instead 
of the correct Pi-q = (M^ + q "^Y^'^E^ — Pj-q in Eq. (|2.8| ); (ii) omit the term l^l,^ in 
Eq. ( p.l2| ); (iii) neglect q^ in the numerators and in the denominators of Eq. ( |2.8| ). The 
first approximation is fatal for Lorentz covariance, the second one for the global variable 
test. 
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B. Modifications of the Riickl formula 



The Riickl formula ([4.1| ) was often used in the calculation of dilepton yield in experi- 
mental and theoretical works. Several modifications of it have been put forward. As will be 
shown later, they brought improvement in a pragmatic sense, but were not able to cure its 
principal drawbacks. 

Gale and Kapusta replaced the photon energy squared ql (cu^ in our notation) in 
the denominator in their Eq. (5) by a symmetrized combination E{E'^ — M^Y^^, where 
E = E+ + E- is the dilepton energy. The quantity (E"^ — M'^Y^'^ represents the dilepton 
momentum and as such must be equal to the photon momentum |q|, which is in turn equal to 
the photon energy go- The replacement Qq E^E"^ —M'^Y^'^ is thus equivalent to multiplying 
the right-hand side of Eq. ( |4.1|) by uo/E. The same modification was used by Haglin, Gale, 
and Emel'yanov in and also in a part of the paper ||2^ by Cleymans, Redlich, and Satz. 

In a subsequent paper Gale and Kapusta introduced a factor which partially cor- 
rected the soft photon approximation for processes with two particles in the final state. It 
accounts for the shrinking of the Lorentz invariant phase space available to them, which 
results from the emission of a dilepton with mass M and center-of-mass-system energy E*. 
The production of dileptons that would violate the energy-momentum conservation is for- 
bidden. The correction factor is given by (we display it in a simpler form assuming equal 
masses m for the final-state hadrons) 



R{s,M,E*) 



s (s2 — 4m^) 
^ S2 (s — 4m2) 



(4.10) 



where s is the invariant energy available for all final-state particles and 



S2 



s + M"^ - 2E*^s . 



(4.11) 



In a paper of Haglin, Gale, and Emel'yanov, the photon energy squared enters 
the denominator of the quantity \^-J\'ih^cd [see their (3.6), (3.7) or (A9)]. It is forced to 
acquire the value of (i?+ -f- EJ)"^ by the second ^-function in their Eq. (3.4). It induces 
a multiplicative factor of uP'jE'^ on the right-hand side of the Riickl formula (|4.1j ). The 



authors also used the correction factor ([4.101). Recently, Haglin and Gale PB| utilized the 



same modification of the Riickl formula to assess the bremsstrahlung contribution to the 
e^e~ invariant mass distribution in proton-proton and proton-neutron collisions at the lab 
kinetic energy of 4.9 GeV. For the final states with more than two hadrons they modified 
the phase-space correction factor accordingly. 



C. Paper by Craigie and Thompson 

After a thorough discussion of the real photon bremsstrahlung, the authors of turned 
to dileptons. If they had really done what they described verbally at the bottom of p. 129, 
they would have obtained immediately a simple and correct formula for the double differ- 
ential cross section in the momenta of leptons [our Eq. (|2.15|) ] with a little different Ansatz 
for the four- vector J (same as was used later in |T^). Unfortunately, they instead wrote a 
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cumbersome and obviously wrong formula (3.1) on p. 130. The incorrectness of the latter 
can be seen, e.g., from the fact that the quantity Tr{pL}, which enters it, depends on the 
momenta of hadrons via the four- vector J [see their Eq. (3.2)]. But there is no integration 
over hadron momenta in (3.1) [only that hidden in 2qQda/d^q, given by Eq. (2.5)]. The left- 
hand side of formula (3.1) in [^ thus depends on hadron momenta, which makes it unusable 
for evaluating the inclusive dilepton cross section. 

The authors probably tried to express the double differential cross section as proportional 
to the cross section for virtual photon production. But, as we have learnt, this is possible only 
for the dilepton cross section in global dilepton variables [compare our ( |2.15D and ( ^.19| )]. 



D. Formula used by Goshaw et al. 



The experimentally observed production of very-low-energy e^e~ pairs in 18 GeV/c 
TT^p collisions was reported and compared to the expectations based on the leading-term 
bremsstrahlung calculations in the paper ||^. The authors used the formula 



da 



a 



[/i2 - (2m 



,2il/2 / 2 2\l/2 



(i/i duj dVt dVt* (27r) 



d^P.---d^P„ 



-{J- J) 



day, 



d^P^ d^P. ■ ■ ■ d^Pr, 



(4.12) 



where /x, cu, and dVL are the dilepton mass, energy, and infinitesimal solid angle, respectively, 
dVL* is the infinitesimal solid angle of positron in the dilepton rest frame, I = P+ — P^, and 

Q^ 



J 



-p 



^(p+ + p_)-i^. 

The charge quantum number of outgoing particles (i = 3, ■ ■ ■ , n) is Qi, of incoming 
ones {—Qi). After changing the notation used in [Q to ours, Eq. ( [4.12|) reads 

,2 ' 



E 



d' 



a^ 



dM^d^q dn+ 327r4 




{i-JgY-Ji 



G 



dan 



where now 



^G 



Pi 



Pi-q 



(4.13) 
= 1,2) 

(4.14) 
(4.15) 



The formula ( 4.14 ) differs from our Eq. ( p.l8 ) by the sign of one of the terms in brackets. 
It is difficult to trace the origin of this discrepancy. The derivation of the formula ( [4. 141 ) 
has not been published, although it was signalized in |^. We suspect that Eq. ([4.14|) is 
not correct, because after integrating it over the positron momentum directions we get the 
formula 



E 



d^a'' 



a 



dW^d^q 67r3 M2 




-Jg) 



(4.16) 



which does not comply with the global variable test, defined by Eqs. ( |1 . 1| ) and ( |1.3| ). In this 
respect, the fact that the authors used additional approximations to get their Eq. ([4.1 3|) 
[compare ( p.8| ) and ( |4.15| )] seems to be of lesser importance. 
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E. Formula of Balek, Pisiitova, and Pisiit 



In paper the formula analogous to our was written for a charged particle 

scattering on a neutral particle (or on a potential). Their formalism satisfies both the 
Lorentz covariance and global variable tests. However, instead of the vector given by 
Eq. ( p.8|) , the following one was chosen: 

0" = ^ ^ . (4.17) 

Pi-q Pa-q 

Here, pa and pi are four-momenta of the charged particle before and after the scattering, 
respectively. In this case, Eq. ( p.8|) becomes 



^ (2pi + qr _ {2pa - qY 

2pi-g + M2 2pa-q-M^ ' ^ ' ^ 

To get the Ansatz ([4.17|) of Balek et al. we have to neglect q^ in the numerators and in 
the denominators above. 

For later convenience we write here the n-final-particle generalization of the Balek et al. 
formula in our notation. The same conventions as used in ( [^.23D apply. 



We will return to the approximation of Balek, Pisutova and Pisut in Appendix 0, in 
connection with the concept of transverse-mass scaling (see below). 



2 



F. Formula of Cleymans, Goloviznin, and Redlich 



The authors of |T3[ used the following classical-electrodynamics motivated expression for 
the energy per unit of momentum radiated in the form of virtual photons with mass M if 
the charged particle changes its velocity from vi to V2 (we switch from their notation to 
ours) : 



d?q 



a 

4^ 



n X V2 



n X vi 



E' — V2-q — Vi-q 



(4.20) 



Our expression for this quantity stems from the relations 



d?q 



E- 



d?q 




d?q 



(4.21) 



M 



where A^'''* is the mean number of virtual photons with mass M per a collision, in conjunction 
with Eqs. ( |2.21| ) and ( |2.2CI| ). Choosing the cross section da^ that allows only the required 
change of velocity, we get 
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with 



V2 + q/(2E2) Vi-q/(2Ei) ^^^3) 



E-V2-q + MV(2E2) Vrq-MV(2Ei) 



Comparing ( 4.2CI|) with ( 4.22 ) and (|4.23| ) we can see that in [|T3|| two additional approxima- 



tions have tacitly been made in contrast to our formalism: (i) The terms proportional to 
q in the numerators and those proportional to in the denominators of Eq. ( [4.23|) have 
been neglected. This is equivalent to the approximation made in [^]. (ii) The second term 
in the brackets of ( [4.22|) has not been considered. This approximation is more dangerous, 



since it makes the quantity (2]20|) non-covariant under Lorentz transformations. 

G. Transverse mass scaling 

In this and next subsections we are going to report about two attempts to relate the 
virtual photon cross section to the experimentally accessible cross section for real photon 
production. 

Farrar and Frautschi W^ and Cobb et al. PP[ used the concept of transverse mass scaling. 



They assumed that the virtual photon cross section does not depend on three variables 
(longitudinal momentum ql, transverse momentum gy, and virtual photon mass M), but 
rather on only two [ql and transverse mass Mt = {ql + M2)1/2]. The condition (|L3D then 
leads to the relation 

where pi, = and pr = (gf^ + M^)^/^. Because both the energies and longitudinal momenta 
of real and virtual photons are now equal, the relation (|4.24]) is covariant under longitudinal 



Lorentz boosts. In our test process —>■ 7r+7r^e+e~, this approach will be very successful. 

The origin and limitations of the transverse-mass scaling from the point of view of the 
leading term bremsstrahlung formalism are discussed in Appendix 

H. Real photon approximation 

Blockus et al. conjectured that the formula 

(j^. d'a^' \ _ fiM,qlqr) ( d^a^ \ 

is valid in the center-of-mass frame (we use ql instead of their x = 2g2/y^). The quantity 
/ is the structure function for the production of a virtual photon of mass M summed over 
photon polarization states. The authors considered two options for their ratio entering the 
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right-hand side of Eq. ( [4.25| ): (i) independent of M (i.e., identically equal to 1), and (ii) 



linearly dependent on M. None of these options seemed to be excluded by their data inte- 
grated over the region of acceptance. In our numerical comparison of various bremsstrahlung 
formulas we will explore the former option 

and will refer to it as the "real photon approximation" . 



V. A MODEL OF TWO-PION RADIATIVE DECAYS OF p° 

As noted earlier, we will check the reliability of the various leading term virtual 
bremsstrahlung formulas by applying them to the strong-interaction decay p° — >■ 7r"'"7r~ 
in order to get estimates of the differential decay width for 

p° ^ TT+TT-e+e". (5.1) 

While a comparison of the bremsstrahlung formulas themselves is instructive, additional 
insight may be gained by also comparing them to results from a formalism that goes beyond 
the leading term approximation. To our knowledge, nobody has investigated the decay ( ^.1|) 
theoretically yet, probably because the chance to detect it experimentally is very meager. 
We present some estimates in Subsection B. 

The real photon counterpart of (|5.1| ), namely p° — 7r+7r^7 has been dealt with by several 
authors pT|j31| , |32[] . In this work we will adopt Singer's approach [^, which agrees nicely 



with the experimental data (see below). We will first recapitulate its main points and then 
generalize it to the massive photon (dilepton) case. 

To simplify notation, we introduce the following ratio of decay widths 

= . (5.2) 

We will call it a branching ratio in spite of a small incorrectness this introduces (FpO^Tr+Tr- 
is smaller than the total width of p° by about 1 %). The quantities B'^* and B^^'^ will have 
analogous meaning. 



A. Decay tt+tt 7 

Strong interaction dynamics enters Singer's calculation through the assumption about 
the p^vr+TT^ vertex in the form 

V'' = f,Upi-p2r , (5.3) 

where pi (^2) is the four-momentum of 7r+ (tt^). For external p's, this is to be contracted 
with the polarization vector ea(Ap). The coupling constant fp-j^T^ can be fixed by utilizing the 
p ^ TTTT decay width. 
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When electromagnetic interactions are switched on by the minimal interaction principle, 
the vertex (|5.3|) generates a contact term 



(5.4) 



which must be considered in addition to the usual vertices of pseudoscalar electrodynamics 
(see Fig. g). It should be contracted with both the p and 7 lines which enter the contact 
vertex together with the two pion lines. Combining the resulting expressions for — > vr+vr" 
and p° 7r^7r^7 decay widths from |]2T]], we obtain the formula 



dui* 



3/2 



2m„ — 2mnU!*] In 



[rrir. 



2uj* 



2uj* 



where lo* is the photon energy in the p° rest frame and 



(5.5) 



2«-^*) 
rup — 2uj* 



(5.6) 



The maximum value of a;* 



is < = {ml-4ml)/{2mp). 
After integrating Eq. ( p. 51) over photon energies greater than 50 MeV we get the branching 



Also, 



ratio of 1.12 %, in nice agreement with the experimental value (0.99 ± 0.16) 
the distribution in photon energies (depicted in Fig. ^ by solid curve) agrees remarkably well 
with experiment [p^ . 



B. Decay —>■ vr+vr e+e 



As has already been stressed, in order to know the differential decay width in global 
dilepton variables of the decay (|5.1|) , it is sufficient to calculate the differential decay width 
for p° — ^ 7r+7r^7* with a massive photon. As was proven in |T^, the evaluation of the latter 
is governed by the same Feynman rules as in the case of a real photon. What changes are 
the kinematical relations which must accommodate the non-vanishing photon mass. As a 
result, we obtain 



dq* 



Aa 



M 



Ami 



3/2 ^* 



E* [m 



2m: 



2E*mp)+M'{E*^ + ^ 



E*q*i 



[mp-2E*J {2mpE*^-M^ 
2 (e*"^ - i'^q*'' 



2m, 



In 



E* + ^q* 
E* - ^q* 



(5.7) 



Here, 
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is the maximum value of E*, the massive photon energy in the p°-rest frame and 



2mp {E^ - E* 



\ml + M^- 2mpE* ' 



(5.9) 



The interested reader may find a few intermediate steps in Appendix Integrating the 
differential branching ratio for p° Tr+vr^e^e^ 



dW^dq* 



'dB^' 
dq* 



(5.10) 



M 



7C ' Tc e ' e 



over the dilepton momenta and masses, we get the value 1.5 x 10 ^ for p° 
and 4.9 x 10""^ for p° ^ vr+vr^/i+^u". A low-dielectron-mass cut of 50 MeV/c^ (100 MeV/c^) 
reduces the branching ratio to 1.1 x 10~^ (4.0 x 10~^). 



VI. BREMSSTRAHLUNG FORMULAS FOR TWO-PION RADIATIVE DECAYS 

OF 

Our goal here is to calculate the differential branching ratio of p° Tr+vr^e+e^ within all 
approaches to dilepton bremsstrahlung we found in the literature. Because some approaches 
||T^ , p!3[] do not provide the distribution in the momenta of electrons, we will concentrate on 
a less general distribution in global dilepton variables. For later use and reference we first 
explore the simpler case of the p° decay to a dipion and photon. 



A. Decay — > vr^vr 7 

In order to calculate the branching ratio in the leading term bremsstrahlung approxima- 
tion, we start with the formula 



duj* Ati"^ 



-4)dQ^* 



(6.1) 



which is a real photon mutation of Eq. ( p.42| ). Using Eqs. ( p.21| ) and ( |2.23| ), we can write 
the integrand in the form 



sin^ a* 



\u* J (1 - v*^cos^a*f 



(6.2) 



where 



1 - 



Am? 



(6.3) 
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is the speed of pions and a* the angle between the photon and vr"*" momenta in the p° rest 
frame. An elementary integration gives us 

dB^ 2a fl + v*^, 1 + v* 
— — = m 

duj* TCUJ* \ 2v* 1 — V* 

A numerical evaluation shows (Fig. ^, dashed curve) that the leading term bremsstrahlung 
formula ( |6.4D exceeds the data, especially at large momenta. This trend is understandable, 
as the formula ignores the energy-momentum-conservation constraints. The branching ratio 
for a;* > 50 MeV is 2.05 %, almost twice as much as the experimental observation. For later 
reference, we also present in Fig. ^ the ratio of the leading-term-bremsstrahlung branching 
ratio ( |6.4| ) to that calculated from Singer's model ( |5.5| ) as a function of the photon momentum 
q* (which is, of course, equal to its energy u*). 




B. Decay — > vr+vr e+e 



Our aim here is to derive formulas for the differential branching ratios in dilepton mo- 
mentum and mass d^B'^*^ /dM'^dq* in the various formalisms. We will apply the same 
assumptions and approximations which the various authors have used when deriving their 
formulas for the differential cross section. For the formalisms that satisfy the global variable 
test ( [|1^,|1^], our formalism), it is sufficient to calculate the branching ratio for the virtual 
photon production 



'dB^ 
dq* 



*2 



M 



a q 



.P] dQr 



[cf. (|2.42|) ] and then multiply it by the function T{M'^), given by Eq. ([T 
the procedure will be less straightforward. 

Applying our formula (p.8|) to the process considered here, we can write 



(6.5) 

In other cases. 



2m~^ 



{E* + M'^/mpY - {v*q* cos6*f 



{E*mp + M^) V * + v*q* cos 6* q * 



Now we insert this into (|2.21|) and calculate the integral which enters the formula (|6 
finish with 



(6.6) 
I). We 



d'^B' 



dM'^dq* 



2a q* 



l + v*'^ - M'^/ml , E* + M^/rrip + v*q* 



2v*{E* + M'^/rUp) E* + M'^/rUp 



V q 



E* 



{v*q*Y + {2E* + M'^/mp)M'^/mp 



(6.7) 



This is the leading term bremsstrahlung formula based on the formalism we developed in 
Sec. fT[ The ratio of ( |0|) to the "exact" formula ( |5.10| ) is presented as a function of dielectron 
momentum q* for several dielectron masses in Figs. ^ through ^ by a solid curve. 

The corresponding formula in the Balek, Pisiitova, and Pisiit approximation is ob- 
tained along the same lines. The only difference is in using their ( |4.17| ), which, translated 



from the scattering case to the p° 



decay, reads 
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2E* 

E*"^ — {y*q* cosO*)"^ 



(6.8) 



instead of our (|6.6|) . The result is 



dM'^dq* 



2a q* 



l + ^;*2 E* + v*q* (l-v*^)E*q* 
m — 



2v* 



E* 



v*q* 



E*^ - {q*v 



*\2 



(6.9) 



It is a good check that for M <^ mp, Eq. ( |6.7|) tends to agree with ( |6.9|) . 

When we combined the Mp-scahng hypothesis ( [4. 24] ) with the leading term 
bremsstrahlung formula for real photons ( p.22| ), we got the same result as from the 
Balek et al. approximation (|6.9|) . This prompted us to investigate the connection between 
the two approaches in more detail in Appendix 0. 

In order to calculate the differential branching ratio in the Cleymans, Goloviznin, and 
Redlich approximation [|1^], we also use (|6.8|). In addition, we have to discard the dot 
product in Eq. (|2.21| ). In this case we get 

^2^e+e- 2a 1 



dM^dq* 



vr E 



E*^ + (q*v*f^ E*+v*q* 
In 



2v*E*q* 



E* — v*q* 



(6.10) 



To find the branching ratios in remaining approximations no additional calculations are 
needed; we need only to combine the proper formulas. 

The Riickl approximation formula ( [4.7|) in terms of branching ratios may be written as 



E 



d'^B^ 



a 



Merging this with 
d^B' 



dW^dq* 
, we get 

a2 1 



27r M2 



M2 ^ dq* ■ 



(6.11) 



dW^dq* 7r2 E*NP 



1- 



4/i2 (l + v 



*2 



2v* 



In 



l + v* 

1 — V* 



(6.12) 



The Gale and Kapusta modification of the Riickl formula consists in multiplying 
( |6.12| ) by the factor q* / E* . To account for the later modification by the same authors |p!T| , 
we must also include their phase-space correction factor Rlrrij^, M, E*), see (|4.10|) . 

To adopt the Haglin, Gale and Emely'anov modification of the Riickl formula, we 



have to multiply ( [TOD by R{ml, M, E*){q* / E*f 



When dealing with the Goshaw et al. 
of (16.51) to the form 



formalism, we first rewrite Eq. ( |4.16| ) by means 



d'^B'^' 



2a 1 



dM'^dq* Svr 



4^2 



dB^' 
dq* 



(6.13) 



M 



The virtual-photon branching ratio on the right-hand side can be taken from the Balek et 
al. formalism, because they use the identical four- vector J [compare ( |4.15| ) with ( [4.17| )]. It 
can easily be read off the Eq. ( |6.9| ); one simply ignores the function T(M^). As a result, we 
get 
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dM^dq* 37r2 (E*M)2 

- v*^) E*q* 




l + v*\ E* + v*q* 
m 



2v* E* — v*q* 



(1 



_ 

The real photon approximation ( [4.26| ) combined with ( p.lOp and (|6.4|) leads to 



(6.14) 



All the formulas we have derived here are normalized to the "exact" formula [Eqs. ( ^.71) 
and ( p.lOl )] and visualized as functions of dielectron momentum at fixed dielectron masses 



in Figs. ^ through | 



VII. COMMENTS AND CONCLUSIONS 



We based our formalism for dilepton production via virtual photon bremsstrahlung on 
the following assumptions: 

1. Only the Feynman diagrams in which a virtual photon line is attached to one of the 
external legs are important. 

2. The charged particles that participate in the process are pseudoscalar mesons or un- 
polarized spin-one-half fermions. 

3. The dilepton four- momentum in the argument of the Dirac 5-function in Eq. ( p. 141) (or 
in similar relations for decays or processes with more than two particles in the initial 
state) may be neglected. 

4. The modifications of electromagnetic interactions of hadrons induced by form factors 
and anomalous magnetic moments (in the case of fermions) are negligible. 

To obtain the leading term approximation, we further neglected terms that are pro- 
portional to the derivatives of the nonradiative matrix element squared. Our notion of the 
leading term approximation is thus based on the proportionality of the bremsstrahlung cross 
section to that of the nonradiative reaction rather than on the order in dilepton momentum 
q. In our leading term formulas, we use the four-vector J that contains also subleading 
terms in q. The numerical comparison with the "exact" formula shows that it is beneficial. 



In Sec. I^, we have pointed out the additional assumptions which are required to reduce 



our formulas to the virtual bremsstrahlung formulas that were derived or suggested previ- 



ously. Of these formulas, only one, namely that of Balek, Pisiitova and Pisiit [|T2[, satisfied 
both the Lorentz covariance and global variable tests. 



The numerical analysis presented in Sec. ^ shows that all leading term virtual 
bremsstrahlung formulas (including ours) overestimate the dilepton production, albeit to 
different degrees. As a matter of fact, also the real photon production is overestimated by 
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the standard leading-term formula (see Figs. ^ and H). But due to the more complex nature 
of virtual bremsstrahlung, the situation with dileptons is more complicated. Besides the ex- 
cess in the high-momentum region, which is similar to that for real photons, some formulas 
also overshoot the yield at small momenta. These formulas relate dilepton production to 
the real photon cross section rather than to the virtual one. Near threshold, the differential 
cross section in virtual photon momentum behaves like g*^ [in both the "exact" formula 
( p.TD and Eq. (|6.7D], whereas for real photons the behavior is q*~^ [see Eqs. ( [5.5D and (|6.4|)]. 

Energy-momentum conservation is not built into the leading term bremsstrahlung for- 
mulas. The fact that they overestimate the dilepton yield at high dilepton momenta is 
therefore not surprising. But we will see that it is only a part of the story. 

Gale and Kapusta fll]] introduced the factor (^4.101) , which accounts for the shrinking 
of the final-state phase space and prevents violation of energy- momentum conservation. 



They, and also the authors of a later work [|J], combined it with the Riickl formula. This 



procedure obscures the conclusions somewhat, because the Riickl formula contains a wrong 
numerical factor. We therefore apply the correction factor ( [4.10| ) to the leading term virtual 
bremsstrahlung formula ( |6.7| ). The results are displayed in Fig. |T0| after being normalized to 
the "exact" formula. Comparison with corresponding non-corrected curves shows that only 
a part of the excess over the "exact" formula has been removed. The remaining excess should 
be ascribed to other sources. It can be a modification of the nonradiative matrix element 
combined with a destructive interference between radiation from the external particles and 
from the contact term. Both these effects are neglected in the leading-term approximation. 

The bremsstrahlung calculations play an important role in assessing the conventional 
sources of photons or dileptons in experiments aimed at revealing anomalous production of 
electromagnetic probes as a sign of new physics phenomena. For example, Haglin and Gale 
have recently shown that bremsstrahlung is the largest source of low-mass dielectrons in 



4.9 GeV pp collisions. A quite good, even if not perfect, agreement with experimental data 
was achieved. The authors used the modified Riickl formula, discussed in [IV B| . But our 
toy example showed that this formula overestimates the dielectron yield by a factor of 2-3. 
This suggests that calculations using a more correct bremsstrahlung formalism have to be 
performed before a definite conclusion about the physics behind the Dilepton Spectrometer 
collaboration data ||3^ is drawn. 

To conclude, we have derived formulas for differential cross sections of dilepton produc- 
tion via virtual photon bremsstrahlung from pseudoscalar mesons and unpolarized fermions 
in the leading and subleading approximations. These formulas satisfy the conditions imposed 
by Lorentz covariance and gauge invariance. From the practical point of view, the leading 
term formulas ( |2.15D , ( [2.1^ ), and ( 2.iy| ) are the most important. They enable us to estimate 
the radiative cross sections on the basis of the corresponding nonradiative cross section. 
Comparing to previously published formulas, our formalism exhibits the best agreement 
with the exact formula when applied to a concrete physical process. 
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APPENDIX A: Mr-SCALING AND VIRTUAL BREMSSTRAHLUNG 

FORMALISM 

Let us describe the multiparticle production in a frame where the velocities of the two 
incident particles are coUinear and define a collision axis. The dot product between the 
dilepton momentum q and the momentum pj of the ith particle can then be written as 

q-Pi = |q|pi,LCOs6' + qr-Pi.T , (Al) 

where 6 is the angle between the dilepton momentum and the projectile velocity. 

If the transverse momenta of the outgoing charged particles are negligibly small (the 
criterion of negligibility depends, of course, on the value of 6), the rightmost term in (|A1|) 
can be omitted. As a consequence, the virtual photon bremsstrahlung formula in the Balek 
et al. approximation ( [4.19| ) can be rearranged into 




^ ^ ~ V^^ l-E-XL|q|cos0 ' ^^^^ 

with Vi^L = pi^L/Ei. Let us introduce a new vector p by 

PL = Ql 



Pt = — ^Jql + M^ (A3) 
qT 



and denote its polar angle by a. We have 



PL |q| a (\A\ 
cosa = - — r = — cost^ . (A4) 

IpI E ^ ' 



Equation 




« f /v-^i Pi 



47r2p2 J y-^ Ei 1 — Vi^i cos a 




(A5) 



The right-hand side is nothing else but the invariant cross section ( |2.22| ) for producing a real 
photon with the momentum p and energy uj = |p| in the same approximation (negligible 



transverse momenta of hadrons). We have thus recovered the Mr-scaling hypothesis ([4. 241) 
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The My-scaling is clearly an approximate phenomenon. To arrive at it, we first made 
the approximations recapitulated in Sec. VII , that led to our basic formulas of Sec. ||. 
Then we added another approximation to get the Balek et al. formula, which is for negli- 
gible transverse momenta equivalent to merging M-r-scaling with leading term formula for 
bremsstrahlung of real photons. In realistic situations, however, the violation coming from 
nonvanishing transverse momenta of hadrons also plays an important role. In our simple 
example-the nonradiative decay —>■ tt+tt" discussed in Sec. [VT] , the outgoing hadrons do 
not have any transverse momenta and the two approaches are fully equivalent. 



APPENDIX B: "EXACT" FORMULA FOR THE p° ^ vr+vr 7* BRANCHING 

RATIO 



The matrix element is given by the FejTiman diagrams depicted in Fig. |^ with a real 
photon replaced by a massive one. We define the invariant variables 



s = (p^ + p^y , 

t' = iP2 + qf - ml 
u' = {pi + qf - ml 



(Bl) 



with pi, P2, and q being the four-momenta of 7r+, tt , and 7*, respectively. The invariants 
( [Bl| ) satisfy the usual relation 

s + t' + u' = ml + . (B2) 

The sum over the and 7* polarizations of the matrix element squared is equal to 

2 



pTTTC 



t'u' 



(B3) 



with the tensor A given by 

A'^" = t'{2p, + qf{pi -P2 + qT + u'{2p2 + qY{p2 - Pi + qT - 2tV(7'^" . (B4) 

Thanks to the relations q^A^"" = and A'^°'Pa = 0, where P is the p^ four-momentum, 
the sums of products of polarizations vectors in (|B^ ) can be replaced by the corresponding 
metric tensors. After a little algebra, we get 



Y: \Mf = lQnafl^ 



Ap,A^* 



2 + "^'^ + ^ 



(B5) 



where 



E,; (mj - 2ml - ^mpE*) + {E*^ + m,/4) 
a2 = \ {ml - 4ml) (4m^ - M') . 



(B6) 
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See also ( ^.81) . To calculate the invariant decay width of the unpolarized p°, namely 



E 



TT 



M 6mp 



d^Pi 



2Ei{2n) 



(B7) 



we need integrals of the type 



fd^Pi fd^p2 1 



(BS) 



which can be computed most simply in the two-pion rest frame after the substitutions 
Q = Pi + P2 and R = pi — p2- We get 



/o = , 
h 



n , E* + ^q* 
m 



rUpq* E* — C,q* 



(B9) 



with ^ given by Eq. (|5.9|) . Putting it all together and using the formula ^1 

P 



3/2 



487rm 



(BIO) 



we eventually get (K 



APPENDIX C: LEADING TERM BREMSSTRAHLUNG AND EVENT 

GENERATORS 

The attitude to the bremsstrahlung in hadronic reactions has been much influenced by 
the fact that the leading term in quantum electrodynamics is equivalent to the corresponding 
classical expression [^]. Similar relation exists also for virtual bremsstrahlung |1T2| . 



As a consequence, also the event generators constructed so far were, according to our 
knowledge, "classical." The configuration of hadrons in momentum space was considered 
a source of photons or dileptons. The momenta of hadrons were assumed untouched by 
the creation of a photon or dilepton. This led to problems with the energy-momentum 
conservation. 

We think that those problems can be cured by a more quantal approach. One should 
consider the transition probability from the initial state to the "complete" final state, con- 
taining both hadrons and a photon (or a dilepton). The energy-momentum conservation is 
firmly enforced. As we show below, the leading term bremsstrahlung approximation means 
that the probability of such a transition can be expressed as a product of two terms. The 
first of them is assumed to be independent of photon or dilepton four- momentum, the second 
one contains the momenta of both hadrons and a photon (a dilepton), but in a simple way. 
What is approximate is the transition probability (rephrased into the cross section, decay 
width, or other observable quantities), not the energy- momentum conservation. 
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To illustrate our point in more detail, let us return back to the basic equations. To 
simplify the discussion, we will consider the reaction 



a + 6-^l + 2 + --- + n + 7. (CI) 
Its cross section is given by (we choose the center-of-mass reference frame) 



^jY.\MH2.)',i,^.n-^P.-,)^, nj^. (C2) 



Due to the presence of the four- dimensional 5-function, only the (3n — 1) momentum com- 
ponents are independent. It is convenient to transform the integration region in (|C2|) into 
a (3n — l)-dimensional unit cube. Several such procedures exist in the literature (see, e.g., 
jH). We thus get 

(C3) 



with 



^3n-l, 



El-^lV(6,6,---,e3n-l) . (C4) 



The function / results from the integration over the four dependent variables and from the 
Jacobian of the substitutions 

Pi = Pi(6,6,---,^3n-i), i = l,2,---,n, (C5) 

q = q(^i,6,---,'^3n-i) • (C6) 



The momenta given by Eqs. (^^) and (|C6|) satisfy the energy-momentum conservation in 
reaction (pi]). The above substitutions are not straightforward, they usually require to solve 
an algebraic equation. We refer the reader to the original literature ||35H38|. 

The master equation of the event generator is equivalent to the evaluating of the cross 
section by a Monte Carlo method: 

where the sums runs over the random points uniformly distributed within the {3n — 1)- 
dimensional unit cube. Each such point generates an "event" -a set of momenta of particles 
in the final state. The quantity ( |C4| ) is a weight that is assigned to each "event". If we 
succeeded in finding such a form of substitutions ( |C5|) and (|C^ ) that the weights of all events 
are same, we would have an ideal event generator. 

We have not made any approximation yet, Eqs. ( [031) and ( P4| ) are exact. Now, we are 
again going to consider only the radiation coming from the external particles and to ignore 
the possible changes in the strong "core" of the Feynman diagram if one of the external legs 
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goes off-shell. The sum of the matrix element squared over the photon polarizations is then 
simply equal to 

Y:\Mf = 47ra\Mof{-Ji) , (C8) 

where Jr is given by Eq. ( p.23|) . To relate |A^oP to observable quantities, we need a set 
of hadron momenta that satisfy the energy-momentum conservation for the nonradiative 
reaction In other words, for the sake of evaluation of the integrand in ( P^) we must 

"spoil" the momenta pi a little. This procedure replaces the assumption that the four-vector 
q in the argument of the 5-function in (|2.14|) may be neglected. We can change the hadron 
momenta in many different ways. Here is one of them. 

We first express C,3n~3, (,3n-2, and ^3„_i in terms of q and remaining ,^'s by inverting 
Eq. ( |C6|) . The momenta of hadrons thus become functions of ^i, ■ ■ ^3n-45 and q. Now 
we can define a set of hadron momenta 

= lim Pi(6, 6, ■ ■ ■ , 6n-4, q) , (C9) 



which satisfy the energy-momentum conservation for reaction We can use them to 

evaluate (|C8| ) by means of the purely hadronic matrix element, which is related to the cross 
section of reaction The weight of the kth event thus becomes 



Ana 



The above procedure can be used if the exclusive cross section is given by an analytic formula. 

In a more realistic situation, the bremsstrahlung cross section is estimated from the 
experimental data on the nonradiative reaction ( |2.1D on event-by-event basis. The Nature 
acts as an "ideal event generator" , which generates exclusive sets of hadron momenta pf^ , 
i = 1, ■ ■ ■ , n. Each set is assigned the same weight (Tq. The distribution of events in hadron 
momenta is governed by the matrix element squared |A^o| and the phase space. The role 
of the latter can again be described using the substitutions 

pf = (^1, 6, ■ ■ ■ , e3n-4) , z = 1, 2, ■ ■ ■ , n (Cll) 



that guarantee the energy-momentum conservation in reaction (|2.1| ). For properly chosen 
substitutions (|C11| ), the ^-dependence of the partly integrated Jacobian /*^°) compensates 
the p-dependence of the matrix element squared in the sense that the product of them is 
constant. The weight of each event thus is, as required, 

^0 = ^^|-Mo|V^°Uei,---,e3n-4) . (C12) 



The momentum distribution of experimental events corresponds to a uniform distribution of 
points within the (3n— 4)-dimensional unit cube. Their coordinates are given by the inverse 
substitution 

6 = 6(pf\pf,---,pi°0 ' ^ = l,---,3n-4. (C13) 
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To generate a photon, we proceed in two steps. Firstly, for each event, we have to "spoil" 
the hadron momenta and add a momentum of photon in such a way that they together cope 
with the energy-momentum constraints for the reaction (PT|). Secondly, we have to find the 
weight of this photonic event. 

In order to obtain the momenta in the radiative event, we supplement the coordinates 
( P13|) by three random numbers ^sn-s, ^3n-2, and ^sn-i- The "spoilt" hadron momenta and 
the photon momentum are now given by Eqs. ( P5| ) and (|U6D . The weight of an event can 
be found by inserting \Mof from Eq. (|CT^ ) to Eq. ( |CTO|) . 



J 5 • • • 5 S3 



w', = A.a .0 i-Jl) l:'^ ■ (CM) 



{P.qK./(o) (^f ),..., 



Of course, the substitutions (|C5[ ) and ( |C11| ) cannot be independent. They must satisfy 



the condition (|C9D . In actual calculation, we do not usually know the ideal substitution 
( |(J11D for momenta in the nonradiative reaction, and work with a substitution that gives 
fluctuating right-hand side of Eq. (|C12| ). The weights of the photon-producing events ( |C14|) 
should be influenced only little, because they contain the ratio f / f^^^ ■ 

The dilepton generators in either global variables or momenta of leptons can be con- 
structed along the same lines with obvious modifications given by the different number of 
variables and different cross sections. 
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FIGURES 




FIG. 1. Matrix element for dilepton production in virtual bremsstrahlung approximation. 
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FIG. 2. Feynman diagrams for vr^vr 7 decay in the approach of Singer |21|. 
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FIG. 3. The differential branching ratio of — > 7r"'"7r 7 as a function of the photon energy in 
the p^ rest frame. Sohd: Singer formula (|5.5|), dashed: leading term bremsstrahlung formula (|6.4D. 



Data 1 22 1 were normalized to the integrated branching ratio given in 
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FIG. 4. The ratio of the differential branching ratio for — > vr+vr 7 in leading term approxi- 
mation to that of Singer as a function of the photon momentum in the p'^ rest frame. 
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FIG. 5. The ratio of the differential branching ratio for p^ 7r"'"7r~e"'"e~ calculated from various 
bremsstrahlung formulas to that of Eq. (|5.7| ) as a function of the dielectron momentum in the p^ 
rest frame at dielectron mass of M = 10 MeV/c^. Solid line: formula ( |6.7| ); 1: Riickl formula; 2: 
formula of Goshaw et al; 3: Gale and Kapusta ||l^ improvement of Riickl formula; 4: Gale and 
Kapusta improvement of Riickl formula; 5: formula of Balek et ai, identical with mT-scaling 
supplemented with the real photon bremsstrahlung formula; 6: formula of Cleymans et al; 7: 
Haglin et al. improvement of Riickl formula; 8: real photon approximation 
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FIG. 7. Same as Fig. | but M = 100 MeV/c^. 
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FIG. 8. Same as Fig. | but M = 150 MeV/c^. 
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FIG. 9. Same as Fig. | but M = 200 MeV/c^. The Riickl formula curve is off the scale. 
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FIG. 10. The role of the phase-space correction factor ( [4. 101 ) in virtual bremsstrahlung at three 
different dilepton masses M: 10 MeV/c^ (upper), 100 MeV/c^ (medium), and 200 MeV/c^ (lower). 
The leading term bremsstrahlung formula ( |6.7| ) with (solid) and without (dashed) the correction 
factor was divided by (|5lo| ). 
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